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We consider nonlinear spectroscopic effects - interaction-enhanced double resonance and spectrum 
instability — that appear in ultracold quantum gases owing to collisional frequency shift of atomic 
transitions and, consequently, due to the dependence of the frequencies on the population of various 
internal states of the particles. Special emphasis is put to two simplest cases, (a) the gas of two- 
level atoms and (b) double resonance in a gas of three-level bosons, in which the probe transition 
frequency remains constant. 
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I. INTRODUCTION 



As is well-known, interaction of a multilevel quan- 
tum system simultaneously with several resonance fields 
is, under certain conditions, accompanied by nonlinear 
spectroscopic phenomena, coherent population trapping 
(CPN) Q} and electromagnetically induced transparency 
(EIT) H, which are caused by the formation of a spe- 
cial "dark" superposition state immune to the resonance 
fields. These effects, however, are still linear in a sense 
that the resonance transition frequencies in the quan- 
tum system are independent of the population of vari- 
ous states. In this work, we consider another, generally 
speaking, a more general class of phenomena induced by 
exactly such a dependence exemplified by collisional or 
contact frequency shift of intra-atomic (e.g., hyperfine) 
transitions in quantum gases owing to the interaction be- 
tween of the gas particles. In our previous work [3|, we 
showed that interaction of a gas of three-level atoms with 
two light fields results, due to the contact shift, in a spe- 
cific kind of double resonance, interaction-enhanced dou- 
ble resonance (INEDOR). In addition, as will be shown 
below, a gas of two-level atoms with a nonzero contact 
shift can exhibit spectrum instability of the resonance 
transition, namely, the dependence of the resonance line 
shape and the final population of the levels on sweep di- 
rection, as well as on the relation between the amplitude 
of the probe field, sweep rate and the magnitude of the 
contact shift. 

The behavior of an arbitrary two-level system is com- 
monly described in terms of effective spin (pseudo-spin) 
s = 1/2. We will also follow this representation, in each 
case attributing spin to a particular pair of quantum 
states coupled by a resonance transition. In this work, we 
disregard spin waves and related effects associated with 
spatial transport of spin polarization, which actually cor- 
responds to zero spin diffusion constant. Below we will 
discuss to what extent this assumption corresponds to 



real physical systems. Spatial inhomogeneity is included 
only in the analysis of a possible INEDOR spectrum line 
shape in a linear gradient of the external field. 



II. NONLINEAR DYNAMICS OF A 
THREE-LEVEL SYSTEM 

In general, the evolution of a quantum system inter- 
acting with resonance fields is described by Liouville-von 
Neuman equation for the components of the spin density 
matrix p [J] 



m^ = [H(t),p]-iC(t) P , 



(2.1) 



where H(t) — Ho + U(t) is the Hamiltonian including the 
unperturbed term Hq and the time-dependent perturba- 
tion U(t) due to the external ac field, C is the Lindblad 
superoperator [f| responsible for dissipation and square 
brackets, as usually, denote quantum-mechanical com- 
mutation. 

To obtain the general evolution equation for the com- 
ponents of the density matrix of a gas of three-level atoms 
interacting with two resonance fields we use the previ- 
ously derived relation for the contact shift of the tran- 
sition |1) — > 1 2) in a spatially homogeneous gas in the 
presence of the third state |3) We will be interested 
in coherent population of the states. In this case, the 
frequency shift Ay(p) = (p) - Wy(0) = 1, 2, 3) of 
the transition between the states \i) and \j) at a certain 
nonzero gas density n = Tr(pjj) vanishes for fermions, 
whereas for bosons in the absence of a Bose— Einstein 
condensate it is 



HA 12 = (pn + p22)<5Ai2 + 

+ (P22 - Pu)AAi2 + 2p 33 (A 

SAi 3 = (pu + p3a)5Ai3 + 

+ (P33 - Pll)AAl3 + 2p 2 2(A 



A13), (2-2) 



A+). (2.3) 
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Here, p^ are the components of the density matrix in 
terms of the eigen wavefunctions of the unperturbed 



Hamiltonian Hn, SXi 
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Figure 1: Scheme of (left) a three- level system and (right) 
influence of the Rabi oscillations between the states jl) and |3) 
on the |1) — |2) transition frequency and intensity. For clarity, 
the Rabi period 27rfi^ 1 of the drive transition is assumed to be 
much longer than the detection time of the |1) — |2) resonance 
line. 



= (ij\X\ij}± is the spin part of the matrix element of 
the interaction intensity A = Airf^a/m, which is com- 
monly used to describe cold collisions, when the par- 
tial amplitudes of scattering with a nonzero angular mo- 
mentum of the relative motion of the colliding particles 
"freeze out" , m is the atomic mass, a is the respective 
s-wave scattering length. The superscript "+" denotes 
that the wavefunction of two colliding bosons is symmet- 
ric with respect to permutation of both their pseudo-spin 
and spatial coordinates. The doubly antisymmetric com- 
ponent, obviously, do not contribute to s-wave scattering. 

Below we restrict ourselves to the typical situation in 
double- resonance experiments, when the state |2) is ini- 
tially unpopulated and its population changes insignifi- 
cantly during the interaction with a weak probe field with 
the frequency w p , whereas the populations of the states 
|1) and 1 3) can vary quite arbitrarily under the action of 
a strong drive field with the frequency uj^- In this case, 
P22 *C pii, P33, and therefore the last term in the right- 
hand side of Eq. (|2.3p can be omitted. Thus, the general 
equation (|2.I[) for the evolution of the components of the 
density matrix of the gas of three-level bosons with a 
ladder (S) level scheme (Fig. [T|) in the absence of spon- 
taneous longitudinal relaxation becomes (cp. [l[) 



.dpi 



dt 



^■pi 3 exp [iAuji 3 (p)t] 



fi, 



P21 exp [iAuii 2 (p)t] — c.c. 



(2.4) 



.dp 



33 



dt 
.dp 3 i 

'~dT 



.9p22 



2 

n 



Pi 3 exp [iAu> 13 (p)t] + c.c, (2.5) 
(P33 - Pxx) exp [iAui 13 (p)t] + 
pp 32 exp [iAu 12 (p)i\ - 1713/031, (2.6) 



%2i exp [iAuJi 2 (p)t] - c.c. 



(2.7) 



dt 



■ dp. 



'32 



dt 



fi P 

~Y {p22 ~ pn)exp[iAuj 12 (p)t} + 

~YP32 exp [iAuJi 3 {p)t] - H12P12, (2.8) 
-■y-pi2 exp [iAu)i 3 (p)t] - 

%-p3i exp [iAuj 12 {p)t] - 1732/932, (2.9) 



where fi p(d) and Aw 12(13) (p) = w p(d) - w 12(13) (p) are the 
Rabi frequency and density-dependent frequency detun- 
ing of the probe (drive) field and 7y are the transverse 
relaxation rates. The light fields are thought to be spa- 
tially homogeneous, which usually holds for hypcrfine 
transitions, whose wavelengths are much larger that the 
geometrical size of the sample. 

If, as already mentioned above, the probe field is much 
weaker than the drive field, fi p <C fid, the second term 
in Eqs. (12. 4p and (12.61) can be neglected. In this case, 
Eqs. (|2.4p - (|2.6[) can be written in a more compact form 
of a usual Bloch equation for the precession of the spin 
polarization ("magnetization") vector M in the Hilbert 
space of the states |1) and |3) (M 2 = (0,0, p 33 — pn), 
Mi = 2(Repi3,Impi 3 ,0)): 



dM 



= m x I n 



AA13 



M z - 7 i 3 Mj 



(2.10) 



where = (fi d , 0,^13(0) -w d + ri 1 SX 13 (p 11 + p 33 ))_. For 
clarity, we separate in the precession frequency Q the 
component that remains nearly constant in a weak probe 
field, when pu + p 33 ~ const. The second term in the 
parenthesis in the right-hand side of Eq. (|2.10l) expresses 
explicitly the dependence of the precession frequency on 
the current value of the magnetization vector M. In the 
subsequent sections, we discuss the direct consequences 
of this, generally speaking, nonlinear precession in two 
limiting cases. 



III. INTERACTION-ENHANCED DOUBLE 
RESONANCE 

As follows from Eq. (|2.2[) . excitation of the transition 
from the state |1) to |3) in a Bose gas induces frequency 
modulation of the transition |1) — |2) associated with the 
Rabi oscillations of the populations of the states |1) and 
1 3) @. If both transitions are excited simultaneously, 
this modulation leads to a new phenomenon, interaction- 
enhanced double resonance. 

Here, for simplicity and clarity, we restrict ourselves 
to the case An = Aj^, when the contact shift of the 
transition |1) — |2) vanishes at p 22 — p 33 = 0. Another 
advantage of this case is that it allows direct compari- 
son with the experiments on electron-nuclear double res- 
onance (ENDOR) in atomic hydrogen [^[l^. To sim- 
plify our consideration even further, we assume that the 
frequency of the drive transition |1) — |3) remains con- 
stant, which requires AA13 to be sufficiently small. An 
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opposite case is considered in Sec. IIVI Despite the above 
assumptions, a general analytic solution of nonlinear sys- 
tem (|2.7[) - ()2.10[) is hardly possible. The approach imple- 
mented below does not feature such generality but pro- 
vides quantitative and physically transparent description 
in the case of interest. 

Physics of the novel kind of double resonance becomes 
clear from Fig. [TJ In the absence of relaxation, the pop- 
ulations of the states |1) and |3) oscillate under con- 
tinuous excitation of the transition |1) — |3) in the ini- 
tially pure state |1) sample at the frequency (see Fig. [T]) 
f2i3 = y/n^+Aujf^. This results in the frequency mod- 
ulation of the probe transition |1) — |2) Q 

hw l2 = hui 12 (0) + 2nAX (j^j 2 sin 2 , (3.1) 

where AA = X 23 — \± 3 . The amplitude of this frequency 
modulation can easily be comparable with or even much 
greater than the linewidth of the probe transition. In 
this case, excitation of the transition |1) — |3) periodi- 
cally drives the probe transition out of resonance. For 
clarity, Fig. [T] corresponds to slow driving in a sense 
that the Rabi period 27rO ( J 1 of the drive transition is 
thought to be much longer than the time needed to de- 
tect the |1) — 1 2) resonance line and the entire |1) — |2) 
spectrum moves periodically fore and back along the fre- 
quency axis and simultaneously changes in amplitude. 
The waveforms shown in Fig. [1] are the schematic "snap- 
shots" of the spectrum at different phases of the Rabi 
cycle. Obviously, electromagnetic absorption at the fixed 
frequency lj p of the probe field also changes periodically 
with a period of the |1) — |3) Rabi oscillations. It should 
be emphasized that these changes are associated with 
changes in both the population of the initial state and the 
transition frequency, in contrast to conventional double 
resonance, which is solely due to a change in the pop- 
ulation of the initial state caused by the drive transi- 
tions. Clearly, such interaction-induced frequency mod- 
ulation can greatly enhance the effect, which therefore 
can be called INteraction-Enhanced DOuble Resonance 
(INEDOR). 

Let us consider in more detail a possible line shape 
of the INEDOR spectrum in a strong drive field when, 
in contrast to the case illustrated in Fig. [TJ simultane- 
ous frequency and amplitude modulation of the |1) — |2) 
absorption line is relatively fast and therefore is inte- 
grated by the detection system. This implies that the 
|1) — 1 3) Rabi frequency is much higher than the rate of 
field or frequency sweep through the |1) — |2) resonance 
or the inverse time constant t -1 of the detection system, 
£l d T > 1. 

Generally, the energies of all three levels and, conse- 
quently, both transition frequencies depend on the exter- 
nal static field. The specific nature of this field is insignif- 
icant but for definiteness we shall consider the magnetic 
field B. Let the field Bq correspond to the exact reso- 
nance |1) — |3) at a certain density p, loi^(p,Bq) = u>d- 



Then the deviation b = B — Bq of the field from this value 
is equivalent to the frequency detuning Awi3(/9, b) = jdb 
and Auj 12 (p, b) — j p b of the drive and probe field, respec- 
tively (here, 7d( P ) is the corresponding gyromagnetic ra- 
tio) . Owing to this Zeeman contribution to the frequency 
f2i3(/9, b) — 7d \J B\ + b 2 of the magnetization precession 
forced by the drive field Bd(t) = B& exp(zwdi), the ampli- 
tude of the oscillating component of iO\ 2 in Eq. (|3.1[) . i.e., 
the amplitude of the frequency modulation of the probe 
transition, is a Lorentzian function of the static field b. 
As a result, uji2(t) oscillates between the Zeeman-only 
zero-density lower bound (dashed line in Fig. [2ji) 

oj 12 {0,B + b) =u 12 (0,B a )+ lp b (3.2) 

and the upper bound, which is the sum of the Zeeman and 
contact-shift contributions (thick solid line in Fig. [2^), 

wia(p, B + b) = Wl2 (0, B ) + lp b + (^^) &2 , 

(3.3) 

at the field-dependent frequency &13Q)) 0|. According 
to Eq. (|3.3[) . the upper bound of the probe frequency is 
generally a nonmonotonic function of the external static 
field. 

The time-average absorption amplitude A(b, uj p ) at the 
probe frequency within the bounds (|3.2[) and (|3. 31) is pro- 
portional to the average population of the initial state 
and the probability density to find the system at given 
values of b and lu p . The probability density, in turn, 
peaks at the lower and upper bound of the probe transi- 
tion frequency, where the system spends more time, since 
in these cases, doJi 2 (b,u> p )/dt = [6|. 

To illustrate the effect of spatial inhomogeneity, let 
us consider in more detail an infinite sample with a spa- 
tially homogeneous density in a linear gradient VB of the 
static field and homogeneous light fields. In this case, 
different parts of the sample simultaneously experience 
all possible field values. The absorption amplitude at a 
given probe frequency lo p is proportional to the integral 
I(uj p ) — j A(b,uj p ) (^) db along the line w p = const 
within the segments AB and CD in Fig. [2^,. The result 
of numerical integration with the parameters of the EN- 
DOR experiments with 2D atomic hydrogen at a density 
of 3 • 10 12 cm~ 2 in the strong polarizing ma gne tic field 
B = 45 kG and the drive field B d = 1 mG 0,113, d, at 
the maximum amplitude AB C = 2riAA(fi,7 p )~ 1 = 89 G 
of the contact shift of the hyperfine transition \b) — \c) 
in field units (here, 7 p and 7d are to a good accuracy 
equal to the gyromagnetic ratio of electron and proton, 
respectively), is shown in Fig. [2Jd a s a function of the 
drive frequency Wd (lower horizontal axis) at constant 
lj p = ^12(0, Bo). Alternatively, the INEDOR spectrum 
can be detected by sweeping w p (upper horizontal axis) 
at constant cjd (in Fig.[5p, = ^13(0, Bo)). In this case, 
the spectrum is inverted on the frequency scale because 
an increase in o>d corresponds to an increase in the res- 
onance value of the static field Wd/7d and, consequently, 
to a positive displacement of the mean-field Lorentzian 
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Figure 2: (a) Field dependence of the |1) — |2) frequency shift (in units of 7 P -Hd) under the CW excitation of the |1) — |3) 
resonance. Horizontal axis is the field detuning (uiis — wa)/7d from the |1) — |3) resonance in units of the excitation field 
Bd- Solid and dashed line are, respectively, the upper bound of the sum of the Zeeman and mean-field contributions and the 
Zeeman contribution alone. Horizontal lines indicate (dash-dotted line) the value of W12 at minimum and (solid line) the probe 
frequency. Vertical dashed line corresponds to the resonance field for the |1) — |3) transition. Inset is the overview of the same 
dependence, (b) |1) — |2) absorption amplitude as a function of the drive frequency (the lower horizontal axis is the detuning, 
in Hertz, from the |1) — |3) resonance) for the monochromatic probe field with the fixed frequency uj p — Ui2(0, Bo). Sharp 
peak corresponds to the minimum of the probe frequency in (a). The spectrum that appears when oj p (upper horizontal axis) 
is swept at fixed ui d = con is inverted with respect to the frequency axis. The parameters correspond to 2D atomic hydrogen 
with a density of 3 • 10 12 cm -2 in the high polarizing field B — 45 kG, except for the sign of the contact shift: Bd = 1 mG, 
AB C = 89 G.' 



peak (|3.3p on the curve W12 (&) in Fig. [5^. The latter is in 
turn equivalent to a decrease in the probe frequency lo p at 
constant Which way of observing the INEDOR spec- 
trum is preferred depends on the details of a particular 
experiment. 

The absorption amplitude as a function of w<j (Fig.03) 
decreases substantially within the drive resonance and 
has a sharp maximum at w p = w m i n , i.e., when the min- 
imum value of Eq. (|3.3I) coincides with the probe fre- 
quency. This explains the dispersion-looking ENDOR 
spectra of 2D atomic hydrogen [9|, [l3| . Qualitatively, the 
hole in the absorption amplitude is because the atoms of 
the otherwise resonant regions of the sample are period- 
ically driven out of the probe resonance and, as a result, 
spend only a small fraction of time in the resonance con- 
ditions. On the other hand, the absorption maximum is 
due to the fact that the drive resonance introduces zero 
gradient of the |1) — |2) transition frequency in a certain 
region of the sample and, therefore, much more atoms 
become resonant. 

The width of the double-resonance curve in the probe- 
frequency units can be estimated as the difference 5u>\2 = 



w m i„ — wi2(0,-Bo) (see Fig. [3^). It is readily shown that 
for a sufficiently high amplitude of the contact shift such 
that |2nAA| > h>y p B d and b > B d 0, 

SoJl3 = 1±Soj 12 ~ ^ 7 d(2AS c i?2) 1 /3 i ( 34) 
7 P ^ 

where AB C — 2nAA(/i7 P )~ 1 is the maximum contact 
shift of the |1) — |2) resonance in field units. Thus, the 
INEDOR linewidth is independent of the static field gra- 
dient. On the other hand, the spectrum intensity is in- 
versely proportional to |V-B|. 

IV. SPECTRUM NONLINEARITY OF A GAS 
OF TWO-LEVEL BOSONS 

In contrast to usual optical Bloch equations, Eq. (I2.10[) 
includes an essentially nonlinear term proportional to 
the contact shift of the resonance frequency. To observe 
this nonlinearity it is sufficient to have just two levels. 
To avoid confusion, we keep levels 1 and 3 omitting for 
brevity the subscripts 13 and p(d). Thus, in this section, 



A A = AA13. In addition, we set the total density con- 
stant, n = pn + p 3 3 = const. The nonlinearity is seen 
most clearly in the spectrum at a low enough sweep rate. 
In particular, if the rate ft -1 AXdM ~ z / dt ~ fi^finAA 
of the variation of the transition frequency due to a 
change in the populations is greater than the frequency 
sweep rate cku/dt, the Rabi oscillations occur on the 
background of a relatively slow change in the magne- 
tization along with the frequency of the light field, so 
that the system stays all the time near the resonance 
Af 2 AA R3 h[uj{p) — ^(0)] — n5X. Meanwhile, the popu- 
lations periodically vary relatively quickly when the res- 
onance conditions are fulfilled, which in turn drives the 
system out of resonance for a certain time until the fre- 
quency of the light field is adjusted to the new value of 
the transition frequency. The picture of such a reen- 
trant resonance repeats until the transition frequency 
stops changing because the limiting value of the mag- 
netization (e.g., M z — n) is reached. Frequency sweep 
in the opposite direction is accompanied by usual be- 
havior of the populations, since the repeated fulfillment 
of the resonance conditions is impossible. On the other 
hand, it is clear that the contact frequency shift affects 
the spectrum considerably if it drives the system out of 
the resonance conditions. This requires that the contact 
shift were at least comparable with the Rabi frequency 
(the bandwidth of the microwave generator is assumed to 
be sufficiently narrow; for example, the spectral width of 
the generator in the experiments with atomic hydrogen 
was Sf/f < 10~ 10 14]) Thus, under the conditions 



ftnlAAl > h 



duj 
~dt 



; n|AA| > Ml 



(4.1) 



there appears a spectrum "hysteresis" (Fig. [3]). If in addi- 
tion n|AA| 3> Ml, a drastic change in the spectrum upon 
reaching a certain critical value of the total gas density 
n c or any other parameter entering the first condition of 
(|4.ip occurs almost abruptly (Fig. [4]) . 

Figure [3] shows the evolution of the population of the 
state 1 3) during the upward and downward (indicated 
by the arrows) linear frequency sweep of the microwave 
field, according to the numeric solution of Eq. (|2.10l) with 
the parameters corresponding to the hyperfine transition 

1 b) — > I a) in three-dimensional atomic hydrogen, namely, 
at the Rabi frequency f2 = 10 s _1 , sweep rate dui/dt = 

2 x 10 3 s~ 2 and the total density n ~ 2 • 10 18 cm" 3 [3], 
AX/h = — 3xl0~ 16 cm 3 /s (see Appendix). The character 
of the spectrum is almost insensitive to the transverse 
relaxation rate (see below). In the calculations, we used 
7 = 0.3 s _1 . The longitudinal relaxation rate in atomic 
hydrogen is quite low and was therefore neglected. The 
total gas density corresponding to each curve is given in 
units of n c m 2.23 x 10 18 cm" 3 . 

The calculated dependence of the final population of 
the state |3) on the total gas density at various sweep 
rates \du/dt\ and the same values of the other parame- 
ters as in Fig. |3J is shown in Fig. 0J When the density 
increases from a subcritical to supercritical value, the fi- 
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Figure 3: Evolution of the population of the state j3) during 
the linear upward (indicated by the right arrow) and down- 
ward frequency sweep of the microwave field according to the 
numeric solution of Eq. (I2.10|) . The total density correspond- 
ing to each curve is marked in units of the critical density 
n c . The values of the parameters correspond to the hyper- 
fine transition b — s- a in three-dimensional atomic hydrogen: 
Q, = 10 s~\ duj/dt = 2 x 10 3 s" 2 , n, 
7 = 0.3 s"\ T = oH, AX/h = -3 x 10~ ie cm 3 /s. 



nal population of the state |3) after sweeping through 
the resonance conditions sharply increases from the low 
value P33 ~ Q 2 \duj / dt\ _1 given by the product of the Rabi 
frequency ft and the duration of the sweep 

through the resonance line, to nearly unity. A further 
increase in density is accompanied by a mere increase in 
the frequency detuning, at which the limiting population 
is reached. As is easily seen, this detuning is exactly the 
maximum possible contact frequency shift nAA. At the 
opposite sweep direction, the final population of the state 
1 3), on the contrary, decreases with an increase in density. 
As one might expect, the final population is independent 
on the sweep direction in the low-density limit, when the 
contact shift vanishes. Thus, the shape of the spectrum 
depends on the rate and direction of the frequency sweep 
of the light field, as well as on the field amplitude. 

As is seen in Fig.[U the critical density n c (determined 
as an abscissa of the steepest slope) does not exactly 
conform to the condition |f2n c AA| oc h\doj/dt\, which fol- 
lows from (|4.ip . and rather increases super linearly with 
\du>/dt\. The origin of this behavior has to be clarifies. 
Note only, that the period 2ir/Q, of the Rabi oscillations 
at a low rate of the microwave frequency sweep (left curve 
in Fig. [4J becomes comparable with the duration of pass- 
ing through the resonance line. Damping oscillations of 
the final population of the state |3) seen at n > n c pre- 
sumably originate from the fact that the system finally 
gets out of the resonance conditions in this or that par- 
ticular phase of effective Rabi oscillations (waves on the 
slanted parts of curves in Fig. [3] at n > n c ). 

Equation (|2.10p holds for a spatially homogeneous sys- 
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Figure 4: Calculated total-density dependence of the final 
population of the state |3) at various frequency sweep rates 
\du)/dt\. Other parameters are the same as in Fig. [3] 



tern. In the case of spatial inhomogeneity, there appears 
magnetization transport owing to exchange and dipole- 
dipole interactions of the atomic pseudo-spins, which 
leads to the emergence of spin waves [1(1 EH ■ The effect 
of the contact shift of the spin-wave spectrum requires 
separate consideration, which lies beyond the scope of 
the present work. 

The effect described in this section can occur not only 
in atomic hydrogen but also in ultracold alkali vapors 
and a number of other systems. In any case, the partic- 
ular character of interaction, which results in the depen- 
dence of the transition frequency on the population of the 
states involved is insignificant. The reason why the spec- 
trum nonlinearity was not observed in the experiments 
on the contact shift in 87 Rb vapor [lj| is that the sec- 
ond condition of (|4.ip was violated. In fact, the s-wave 
scattering lengths of 87 Rb in various hyperfine states are 
such that the maximum differential contact shift at the 
density n ~ 10 13 cm -3 was as small as nAX/h ~ 0.2 Hz 
(Auj/uj ~ 3 • 10~ n ), whereas the Rabi frequency of the 
two-photon transition in the pulsed microwave/RF field 
was about 2.5 kHz. 

The dependence of the resonance field on the sample 
magnetization knowingly leads to a so-called ferromag- 
netic instability [Hj], e.g., of the FMR spectrum of ferro- 
magnetic films, ESR in atomic hydrogen J17| and NMR 
in 3 He In contrast to this type of instability, whose 
condition is determined the transverse relaxation rate, 
the effect considered in this section does not depend di- 
rectly on the transverse relaxation, as the contact shift 



in a two-level system is independent of the mutual coher- 
ence of the single-particle states [19( • 
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cussions and providing us with the data of atomic hydro- 
gen experiments the University of Turku, Finland. This 
work was supported by the Human Capital Foundation, 
contract no. 211). 



Appendix A: Contact Shift in Atomic Hydrogen 

The magnitude of the contact shift of the hyperfine 
transition \b) -f-> \a) in atomic hydrogen can be found as 
follows. The diatomic states in the basis \S, ms', I, mi) of 
the total electron and nuclear spins of the pair of atoms 
take the form 



\bb) 
\ab)+ 

\aa) 



|i,-i; i, 



(Al) 



— Qab) + \ba)) 



cos0|l,-l;l,O) -sin(9|l,0;l,-l), (A2) 
cos 2 0|l,-l;l,l) 1 — 2 ' 
sin 20 , 



|1,0;1,0) 



sin" t>\l, 1; 1, -1) 
sin 20 , 



■|0,0;0,0) (A3) 



where tan(20) = A[{^ e + r y p )hB}^ 1 , 7 e (7 P ) is the gyro- 
magnetic ratio of electron (proton), A/h = 1420 MHz 
is the hyperfine constant of hydrogen). That is, the 
states 1 66) and |a6)+ are pure electronic and nuclear 
triplets irrespective of the value of magnetic field. Con- 
sequently, X^ b = Xbb = a t , where a t is the triplet 
s-wave scattering length. Thus, according to Eq. (12.31) . 
the contact shift of the transition |6) — \a) vanishes at 
Paa = in an arbitrary field. On the other had, the 
state \aa) contains the singlet component (the last term 
in the hight-hand side of Eq. (|A3| ). Consequently, 
AA ab = X bb + X aa - 2A+ = ^(a s - a t )sin 2 20 ^ 0, 
a t - a s = 30(5) pm [ljj. In the field B = 4.5 T, 
AX a b/h = —3 x 10~ 16 cm 3 /s, and the frequency shift 
at the density p aa = 2 ■ 10 18 cm -3 is about —100 Hz, 
which is two orders of magnitude greater than the Rabi 
frequency. 

Replacing the electron spin by the nuclear spin all 
the aforesaid is automatically generalized to the tran- 
sitions | b) O |c), | a) O \d) and \c) ■(-> \d), wherefrom 
AA 6c = AA ad = AA cd = ^(a s - a t )sin 2 20. However, 
the populations of the states with opposite electron spins 
cannot be simultaneously high owing to a fast recombina- 
tion of such pairs; thus, the sift of the transitions |6) ■(-> \c) 
and | a) O \d) is hardly detectable. 
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